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Abstract. This work confirms the stability of a class of domain wall lattice models that can 
produce accelerated cosmo logical expansion, with pressure to density ratio w = —1/3 at early 
times, and with w = —2/3 at late times when the lattice scale becomes large compared to the 
wall thickness. For walls of tension Ti, the relevant X type junctions could be unstable (for 
a sufficiently acute intersection angle a) against separation into a pair of Y type junctions 
joined by a compound wall, only if the tension Tj of the latter were less than 2T\ (and for an 
approximately right-angled intersection if it were less that y/2T{) which can not occur in the 
class considered here. In an extensive category of multicomponent scalar field models of forced 
harmonic (linear or non-linear) type it is shown how the relevant tension - which is the same 
as the surface energy density U of the wall - can be calculated as the minimum (geodesic) 
distance between the relevant vacuum states as measured on the space of field values $ l using 
a positive definite (Riemannian) energy metric dll 2 = Gij d$* d& that is obtained from the 
usual kinetic metric (which is flat for a model with ordinary linear kinetic part) by application 
of a conformal factor proportional to the relevant potential function V. For suitably periodic 
potential functions there will be corresponding periodic configurations - with parallel walls 
characterised by incrementation of a winding number - in which the condition for stability of 
large scale bunching modes is shown to be satisfied automatically. It is suggested that such 
a configuration - with a lattice lengthscale comparable to intergalactic separation distances - 
might have been produced by a late stage of cosmological inflation. 
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1. Introduction 

This article is concerned with the question raised by Bucher and Spergel [TJ [2] of whether 
some kind of cosmic domain wall lattice might account for the observed acceleration of the 
expansion of the universe, not to mention subtle deviations from isotropy[3J H|- More specifi- 
cally, the question is whether the tendency of such lattices to evolve - typically according to 
a scaling law [5] - towards a uniform vacuum solution can be "frustrated" in particular sce- 
narios in which the lattice is ultimately preserved in a configuration that is effectively frozen 
with respect to comoving coordinates. A first rerequisite that must be satisfied by proposed 
candidate models 0, H E] for such a frozen state is of course that of stability, but even in cases 
for which this is satisfied there remains, as a further necessary condition for viability, the 
more difficult problem of attainability of the configuration in question from plausible initial 
conditions. 

The first aim of the present work is to show how the existence of absolutely stable frozen 
lattice configurations will be an automatic consequence of the existence of conserved topo- 
logical winding numbers for an extensive class of field models with multiply connected field 
configuration spaces. Although such models would seem to be capable of providing good 
agreement with what we see now, this depends on the provision of initial conditions for which 
- contrary to what might be expected from naive symmetry considerations - the relevant 
winding numbers are endowed with non zero values. That nature might conceivably provide 
such symmetry violation - perhaps for anthropic reasons - is shown by the notorious example 
of cosmic baryon number asymmetry. The need to invoke such a priori symmetry violation 
does however diminish the attractivity of such scenarios. 

In the (not yet fully satisfied) hope of getting round this drawback, the following work will 
be mainly concerned with the extension of the class of field models under consideration to a 
more general class in which, when sufficient energy is available, the fields will have access to 
a configuration space that is simply connected, but in which the energetically attainable field 
values will otherwise be effectively confined to a neighbourhood that is postulated to have the 
non-trivial multiply connected kind of topology considered in the preceding paragraph. In sce- 
narios based on such models, cooling from a thermally excited state can be expected to leave 
field configurations with winding numbers that will be effectively conserved by substantial 
energy barriers, and that may be left with non-vanishing values - thus guaranteeing "frustra- 
tion" - over causally connected volumes that in an inflationary scenario could be larger that 
the visible universe. However an (in the present context) undesired biproduct of this process 
would be the formation of string type defects on which the domain walls resulting from the 
winding would terminate. It is easy to conceive ways in which such undesired strings might 
be inflated away (as in the usual solution of the monopole problem) but would seem that to 
do this without also inflating away the desired wall lattice would again require recourse to fine 
tuning (and perhaps invocation of the anthropic principle) such as we were trying to avoid. 
The ultimate message of this article is that this approach does not seem very promising at the 
present stage, but that it should not yet be definitively excluded from consideration. 

It is to be remarked that a more decisively negative opinion, namely that viable scenarios 
for a "frustrated" lattice cannot exist, has been vehemently advocated in recent articles by 
Avelino et al. [9j [10]. (Since they went so far as to claim that the viability of my proposed 
pentavac example [6] was "easily" ruled out by their numerical work, despite the fact that the 
latter was confined to a limited - apparently inappropriate - part of the relevant parameter 
space, the present article will therefore give particular attention to the provision of analytic 
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reasonning showing that this pentavac model does indeed provide lattices with the required 
stability properties provided the relevent parameters are chosen within the appropriate range 
as evaluated in an appendix.) In more recent writing [IT] these authors have conceded that our 
examples 0, [7J, E] do show how "one can build (purely by hand) special lattices that would be 
locally stable against small perturbations". They nevertheless maintain their no frustration 
conjecture to the effect that "no such configurations are expected to ever emerge from any 
realistic cosmo logical phase transition". I would suggest that a more justifiable conclusion 
(from their own and other concordant work on the tendency towards scaling behaviour [5]) 
would be obtainable by substituting the qualification "easily numerically simulable" in place 
of their adjective "realistic", but I agree that the "frustrated" examples in this and the pre- 
ceeding work 0, [JJ [8] can quite fairly be criticised as artificially contrived. What should not 
be overlooked however is that one can also criticise wild flowers as having been artificially con- 
trived to attract bees: the point is that artificially contrived results are sometimes obtainable 
by Darwinian or anthropic selection mechanisms in a manner that is undeniably natural and 
"realistic" , albeit far beyond the scope of easy numerical simulation. 

Before proceding to what is new, this article will start by recapitulating some noteworthy 
conclusions from the preceding work [6j [7J [8] in which attention was drawn to the important 
qualtative distinction between the X type of (crossover) junction that is more favorable for 
stability of the lattice, and the Y type of junction, in which there is no freedom of adjustment 
in the equilibrium angle of intersection (which must be ir/3 if the wall tensions are all equal). 
For an X type crossover the most symmetric possibility is a right angle intersection, but 
for opposing pairs with equal tension equilibrium will still possible when there is a positive 
deviation 5 so that the walls meet at an acute angle 



It has however been emphasized [9] that in order to contribute to a stable lattice such an X 
type equilibrium must be stable against decomposition into a pair of Y type junctions (see 
Fig. [1]), a requirement that was not explicitly checked in the particular toy field models I 
suggested [6] as examples in my original discussion of this subject, and that has been called 
into doubt [5] in the particular case of what will be referred to here as the pentavac doublet 
model. 

In order to address the lattice stability question, a preliminary task of the present work will 
to provide a simple general criterion for the stability of such an X type equilibrium junction. 
A strategy for testing this criterion will then be provided for an extensive class of forced 
harmonic field models, including the pentavac doublet model [6], for which it is confirmed 
that (contrary to the doubts that have been expressed [9j[lQ]) the stability condition is indeed 
satisfied. The main part of this work is concerned with the related but more delicate issue of 
stability against bunching of parallel walls, which will be dealt with in the same framework. 

Although - like many other possibilities such as the monovac triplet model developed at 
the end of this article - the pentavac doublet model can provide a regular lattice that is 
stable, it should be mentioned that it is nevertheless unsatisfactory from the point of view 
of the purpose that motivated its introduction, which was to provide a random lattice of the 
kind to which my (still unproven) five colour conjecture was concerned[6]. The special feature 
of the pentavac doublet model (see Fig. [4] at the end) is the admission of simple domain 
walls between any of the 10 pairs that can be chosen from its 5 distinct but equivalent vacua, 
and the admission of X type crossovers involving any of the 5 possible combinations of 4 
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distinct vacuua. However, for each such combination, the pentavac model allows only one of 
the 3 mathematically conceivable ways of choosing the dagonally opposing pairs, whereas all 
of them would be needed for a random solution of the five colour problem. 

This limitation on all the cases investigated so far - namely that they provide lattices 
that can be stable only when sufficiently regular, not random, and hence that they seem to 
depend on prerequisite "tuning" of the universe - has considerably reduced the attractivity 
of such domain wall models in comparison, for example, with the simple hypothesis of an 
appropriately "tuned" cosmological constant. 

Despite this limitation, wall scenarios of the periodic lattice type considered here remain 
viable in principle, as a possibility that should perhaps be taken more seriously, particularly 
[31 H] if the evidence for cosmological anisotropy [12] is confirmed. The systems proposed for 
investigation here are of a kind that would arise from multiplet generalisations of the much 
studied Peccei-Quinn singlet model, which produces axionic walls whose potentially important 
cosmological consequences have been considered by Khlopov and coauthors [131 dH US]. In 
systems of this kind, the collective stability of the walls depends on their feature of having 
topological winding numbers, which need to add up coherently with the same sign. The 
problem with this is that in a random system (such as would be expected from a Kibble type 
symmetry breaking mechanism [T6l [T71 [T8] ) one would expect to obtain roughly equal numbers 
of positive and negative windings (which in the long run would undergo mutual destruction, 
leaving a residue of string anomalies). 

In comparison with the (unavoidable) problem of accounting for the observed (but still 
mysterious [19]) preponderance of ordinary matter over antimatter, the analogous, but numer- 
ically less extreme, problem of accounting for the required preponderance (on a sufficiently 
large scale) of positive over negative walls appears to be less intractable. The reasonning 
developped below suggests that it may be soluble on the basis of statistical fluctuations in the 
framework of suitable inflationary models - albeit with the help of special parameter tuning 
such was already inherent in such models. The basic idea is that - assuming the strings that 
might have been formed in a very early high (e.g. GUT) energy transition were subsequently 
inflated away - reheating could have produced walls associated with more moderate (e.g. elec- 
troweak) energy scales, that (due to statistical fluctuations) need not have exactly cancelled 
themselves out in a localised volume, but could have left a residue with coherent winding on 
a sufficiently large (cosmologically significant) scale. 

2. Simple and compound wall configurations in models with X type junctions. 

The recently raised issue [9] of possible instability of X type junctions arises in bosonic field 
models of the kind I suggested as candidates for providing domain wall domain wall lattices, in 
which the essential feature - as illustrated in Fig. H]- is the existence in the space of classical 
field values of neighbouring subsets of four equivalent discrete vacuua - meaning minima of 
the potential energy V - separated by four ridges of higher energy that meet in a cross shaped 
configuration at a peak of even higher energy. Using the letters A, B, C, D, to label the energy 
minima in cyclic order, as one goes round the peak of V in configuration space, entails the 
concomitant notation AB, BC, CD, DA to label the ridges that separate them. 

For each such ridge in the space of field values, there will be a corresponding flat domain 
wall equilibrium state in ordinary spacetime: for example AB will designate the - energy 
minimising - domain wall state specifiable as a function of a single cartesian cooordinate, 
x, say by the condition that it tends to the field configuration A as x — > — oo and to field 
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Figure 1: Transition from X junction (deviating from orthogonality by angle 5) to neighbouring 
pair of Y junctions (with complementary intersection angle a) for walls between vacuum 
domains labelled A, B, C, D, using heavy line to indicate compound wall. 

configuration B as x — > +00, with maximum energy density at x — 0. Since there is nothing 
to break the Lorentz invariance in directions parallel to the x = plane, such a wall will be a 
Dirac type brane, with isotropic tension, T\ equal to its surface energy density. 

As well as such simple branes associated with the field space energy ridges AB, BC, CD, 
DA, there may also be compound wall configurations separating non adjacent vacuum pairs, 
for which the possible combinations are AC and BD. In order to be stable against splitting into 
the relevant pair of simple walls (for example in order for the compond wall AC to be stable 
against splitting into the separate simple walls AB and BC) the field model should evidently 
be such that the energy density (and tension) Tj say of the compound wall configuration - if 
it exists - is less than twice that of a single wall: 



3. Criterion for stability of an X type junction 

In a model of the kind considered in the previous section, an ordinary X type junction 
between four vacuum domains A,B,C,D, consists of a string like locus of intersection through 
which a simple wall of type DC continues as a simple wall of type AB, while a simple wall of 
type AD continues as a simple wall of type BC. Such an X type junction will always be stable 
if the model does not admit any compound wall configurations satisfying (J2j). 

The possibility [9] of instability with respect to decomposition into a pair of Y type junc- 
tions will however arise if the model is such as to admit compound walls satisfying the condition 
(j2J). If the simple walls AB and BC meet at an acute angle a , the decomposition in question 
would create a Y junction joining them to a compound wall of type AC, while - as shown in 
Fig. [T]- at the other end of the double wall there would be another Y type junction on which 
the simple walls AD and DC would meet at the same acute angle a. The original X type 
junction will be unstable if and only if the double wall of type AC is too weak to prevent the 
separation of the two Y junctions from increasing, that is to say if and only if 



(2) 



T n < 2Ticos{a/2} 



(3) 
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a condition that is evidently stronger than (j2J) This is expressible the other way round as the 
condition that the X junction will be stable if and only if the acute angle a is such that 

cosW2}<^. (4) 

Since the range of possible values of the intersection angle is given by0<a<7r/2, which 
is equivalent to 1 > cos {a/ 2} > l/y/2, it follows firstly that, in the regime for which 

Tn>2Ti, (5) 

(so that the the criterion (j2J) for local stability of an individual compound wall will fail) the 
stability condition (T4J) will always be satisfied. Secondly, in the regime for which 

2T l >T JL >V2T I , (6) 

there will be a critical crossing angle, given by 

rp 2 

cos a c = sin 5 C = — ^ _ 1 > ( 7 ) 

such that the X junction stability condition (J3J) will hold if and only if the actual crossing 
angle satisfies the criterion expressible in the equivalent forms 

a > a c , 8 < 5 C , (8) 

in which when (JHJ) holds, it will be possible to construct a stable periodic X-junction 

lattice of the kind illustrated in Figure [21 Alternatively, if the crossing angle were too acute for 
(jBJ) to hold, the same global boundary conditions could be satisfied by a periodic Y junction 
lattice of the kind illustrated in Figure EJ but - as in previously discussed examples [SI El IS] 
- the stability of such a configuration with respect to local perturbations would be marginal. 
Finally, if 

TzKyfiTi. (9) 

it is evident that the condition (j4j) for X junction stability will never be able to hold at all 
for any value of the crossing angle. 

4. Scalar field models of forced harmonic type 

The toy bosonic field models we have been considering belong to a rather extensive category 
that is describable as being of forced harmonic type. This means that the independent scalar 
field components, $ 4 say, are to be considered as coordinates on a manifold characterised by 
a flat or curved Riemannian (positive definite) metric 

dA 2 = Gij d$* d$ J , (10) 

and by a scalar forcing potential V, in terms of which the relevant Lagrangian density in 
ordinary spacetime, with pseudo Riemannian metric dfi 2 = g^ u dx^ 1 dx u takes the standard 
form 

c = -±G ij (y'&)v v &-v. (ii) 
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Figure 2: Locally stable periodic lattice satifying stability condition or (jSj) with subcritical 
deformation angle S for X junctions between walls separating vacuum domains that might be 
identical, as in monovac winding models, but that might also for example be of 4 types as 
indicated by the labels A, B, C, D, or even of 5 types as indicated by the coloring. 




Figure 3: Periodic lattice for same kinds of model with same global (phase winding) boundary 
conditions and deformation angle 5 as in Figure [2] but in regime (EI) with smaller critical value 
5 C < 5 so that junctions will be of Y type, with critical intersection angle a c = it/ 2 — 8 C . This 
lowers energy density by factor cos{(<5 — S c )/2} , giving configuration whose stability against 
local perturbations is marginal. 
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Generalised non-linear sigma models (of the kind whose properties were used by Bunting [20] 
for his proof of black hole uniqueness) are included as the subcategory for which the potential 
function V is set to zero. Scalar field models of the ordinary linear and non-linear kinds that 
are more commonly considered [51] in physics are included as the subcategory for which the 
forcing potential is a variable function, U{$}, but for which the kinetic metric ( flQj) on the 
space of field values is flat, so in such a kinetically linear model it will be possible to scale the 
field variables in such a way that the metric components will be given simply by the Kronecker 
unit matrix as Gij = Sij. 

In a flat static domain wall state of the kind considered in the preceding section, for which 
the fields depend only on the single cartesian space coordinate x, the surface energy tension 
T will be the same as the integral along the x axis of the energy density, as given in the limit 
t — > oo by the prescription 

T = U{oo}, (12) 

where 

ft/2 ( d^d^' \ 

u ^=L{^^ +v r' (13) 

in which (to get a finite result) it is to be understood that the potential V has been adjusted 
if necessary, by the addition of a constant which will have no effect on the field equations 
obtained from fill I) , in such a way as to arrange for its minimum (vacuum) value to be zero, 
V = 0. 

Stable wall equilibrium states are characterised by the condition that the parametrised 
field manifold trajectory (as specified by the functions $*{x}) with endpoint at the relevant 
pair of vacuum position (where ^{$} is minimised) should be such that the energy integral 
(|12[) should be a minimum. 

Using a prime to denote derivation with respect to the parameter A measuring the distance 
in the field manifold, so that in particular, according to (flOl) . we shall have 

G ij & , &' = 1 , (14) 

it can be seen that the integrand in f[T3"|) can be rewritten as 

dU = U'd\, (15) 



with 



in terms of the rate of variation 



U '=h +Vx '< (16) 

(17 » 



of the space coordinate x with respect to the field manifold distance given by ffTO]) . For 
given field values at the extremities of some finite range — 1/1 < x < £/2, a corresponding 
equilibrium configuration will be characterised by minimisation of the energy integral (fT2]) 
with the integrand given by 

U{£} = Ju'dX, (18) 

subject to the associated constraint 

f x'd\ = £. (19) 
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In terms of an appropriately adjusted value of a Lagrange multiplier P± say, this condition 
will be equivalent to unrestricted minimisation of the corresponding enthalpy combination 



H{£} = U{£} + P ± I (20) 

which will be given by 

H = J Q^ + (V + P±)x') dX. (21) 

By considering the way this depends on x' for a fixed dependence of the fields on the parameter 
A, it can immediately be seen to be necessary for minimisation that the variation rate x' should 
satisfy the relation 

*' 2 = WT^ry (22) 

in which P± is interpretable as a constant of integration whose value can be seen to be that 
of the pressure in the direction orthogonal to the plane of the wall. 

The condition (T22]) can be used to eliminate the involvement of the ordinary space coordi- 
nate x, and to express the quantity that (for fixed field values at the endpoints of integration) 
has to be minimised in the form 

H = J ,J2(V + P ± ) dX, (23) 

in which only field manifold variables are involved, and in which the fixed parameters £ and 
P± are not independent but according to (1221) must be related by the consistency condition 

(24) 



2(V + P ± ) 

It can be seen from this that when the orthogonal pressure is varied the corresponding 
variations of U and I will be related by the formula 

whereby the interpretation of P± as the relevant orthogonal pressure is made obvious. 

In cases such as those of the compact, topologically non-trivial, field manifolds to be 
considered in the next two sections, a solution for a finite value of £ may be extensible over 
the complete range of x as a "lasagne" type periodic configuration of the kind whose (linear 
or non-linear) superposition will contitute the kind of cosmological lattice whose investigation 
is the ultimate motivation for this work. More specifically, the slab thickness £ will then 
be identifable with the wavelength of the periodicity, so that the corresponding longitudinal 
wavenumber density will be 

provided the endponts of the integration occur at successive maxima or successive minima of 
the potential as function of x. Thus in particular, for P± > 0, the distance £ will be the same 
as the wavelenth if the endpoints occur at successive vacuum states. 

In such an effectively one dimensional lasagne type configuration, the condition for stability 
of the large scale averaged system against development of a bunching mode is the positivity 
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the square v]_ of the orthogonal propagation speed, of perturbations of the mean density U u± 
- on length scales large compared with the wavelength £ - as given by 

2 _ dP ± _ £ dH 
V± ~d(Uu ± )~ Hd£- 

That the requirement v]_ > will indeed be satisfied is shown by the expression 

- 2 H dl d£ f fn ,„ . „ xX -3/2 



*?= -iPi=J^ v+p ^ dA ' (28) 

It can be seen (using the Schwarz inequality) that the velocity given by (1251) can never 
exceed unity, meaning the speed of light. It will however approach this upper bound when P± 
becomes very large, so that the kinetic energy becomes much greater that the potential energy 
(whose particular form will then be irrelevant) which is what happens in the short wavelength 
limit i — > 0, for which one will have 

U~p x £, P ± ~I(A/£) 2 , wi ~ 1, (29) 

where A is the integrated distance along the field space trajectory from one vacuum state to 
the next, as measured with respect to the kinetic metric (flOl) 

Our main concern here is not with the short wavelength limit, but on the contrary with 
isolated wall configurations, which are obtained when the range of integration is unlimited, 
C — > oo. This large separation condition requires that the integral (124j) should be divergent 
at the extremities of the trajectory in field space, namely the vacuum states characterising 
the domains in question, where the potential V reaches its minimum. Since, for applicability 
of the tension formula (1T2|) . it is to be understood that (by subtraction, if necessary, of any 
contribution from a cosmological constant that may be present) this minimum value of V has 
been adjusted to zero, it can be seen that to obtain an isolated wall configuration the relevant 
constant of integration must also be taken to be zero, 

P± = 0. (30) 

According to (1T6|) this simply gives 

U' = VW, (31) 
and allows one to take the limit £ — > oo in ( |23l) to obtain an expression of the simple form 

[/{oo} = / VWdX, (32) 



for the surface energy to be minimised in order to obtain the equilibrium configuration of the 
isolated wall. 

It can be seen from the formula (13~T1) that this required surface energy density function 
U is interpretable as a generalisation to the multiscalar case of what has been referred to in 
the context of a singlet field [22] as a superpotential, and that it is specifiable as the measure 
of the relevant distance in the field manifold, as evaluated, not with respect to the original 
kinetic metric Cry, but with respect to a conformally modified field energy metric that is given 
by 

dU 2 = Gij d& d$ j , (33) 
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with 

= 2VGij . (34) 

The geodesic distance (between the relevant vacua) obtained by minimising the integral £/{oo} 
given by (1321) will thus be directly identifiable with the required wall tension T as given by 
(1121) . The corresponding distribution of the fields as a function of the orthogonal distance x 
is given - according to (1221) - by a metric that is related to the kinetic metric by a conformal 
factor that is exactly the inverse of the one that gives the energy metric: 

dX= 2V=Wf (35) 



5. Analytically integrable cases of pentavac and monovac doublet models 

As an illustration of the kind of model that can provide a lattice with stable X-type 
junctions, I proposed as an interesting prototype example the particular case of what may be 
concisely referred to as the pentavac model pj. This is a model with five equivalent vacuum 
states in a toroidal field space, in which the relevant field variables are a pair of phase variables 
4> and if> with period 2tt with flat kinetic metric given in terms of some fixed mass scale, r/ say, 
by 

dA 2 = r] 2 (de 2 + d X 2 ) = 5r] 2 (d(j) 2 + dip 2 ) , (36) 
and potential V given in terms of some fixed maximum value V* by 

V= -^(cos0 + cosx + 2). (37) 

in which 

# = 20 + ^, X = 2^-0. (38) 

A model of this topologically non-trivial kind is obtainable j6] as a low energy limit e — > 
from a topologically simple extended model of with broken U(l)x U(l) symmetry, of the kind 
is decribed in the appendix. 

According to (I3~4"|) the energy metric obtained from flHBI) and fl3T|) on the toroidal space of 
phase variables and if> will be given by the formula 

dU 2 = 2 v 2 V(d6 2 + d X 2 ) , (39) 

in which it is useful to rewrite the formula (I3TI) for the relevant potential as 

V = V* (icos 2 {#/2} + |cos 2 { X /2}) . (40) 

For such a potential, the corresponding geodesic Hamilton Jacobi equation, namely 

1 , (41) 



'dSS\ 2 fdSSV' 



2rj 2 V \ I 89 J \dx 
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turns out to have the convenient property of being separable with respect to the variables 9 
and x ■ By taking the Jacobi action variable SS to be a sum, 



SS = SS e + SS X , (42) 
of single variable functions, the equation (}4"Ij) can be seen to be reducible to the form 



A C O \ ^ //ICC ^ 

(AO OQ \ o Tr 9 r . o Tr 9r / QjJx 



- V ^ cos^/2} = cos^{x/2} - j , (43) 

in which the terms on the left depend only on 9 while those on the right depend only on x, 
which means that both sides must be equal to a constant of integration, K, say. This means 
that the corresponding generalised momentum variables, which are specifiable according to 
(PD as 

-Pe = 2rfV^j , V x = 2rfV^ , (44) 
will satisfy equations of the form 

V e 2 = r] 2 V* cos 2 {#/2} + K , V 2 = tfV* cos 2 {x/2} - K . (45) 



It is to be recalled that the only trajectories to which this derivation is applicable are 
those that begin and end on vacuum states, namely the field values for which cos 2 {0/2} = 
cos 2 {x/2} = . In order for the momentum values given by fj4"5"]) to remain real at these 
end points it is evident that neither K, nor —K can be negative, which is only possible if the 
constant itself vanishes, 

K = . (46) 

It can thereby be concluded that the equations of motion for the relevant geodesies are given 
simply by 

4^ 2 (^) 2 = K cos 2 {tf/2} , A7] 2 V 2 f~^V = V * cos 2 {x/2} . (47) 

It can be seen from this that, independently of parametrisation, the trajectory in phase 
space will be obtainable by integrating 

d " ■ d * (48) 



cos{#/2} cos{x/2} ' 

and that the corresponding expression for the energy variation (in the "right" direction) will 
be given by 

dU = 2r] yjv*(caa{9/2}dd ± cos{x/2}dx) . (49) 

An example of a simple wall trajectory (such as the straight line AB in Fig. H]) is obtainable 
by holding \ fixed at a value such that cos{x/2} = (so that the ratio on the right of (|4*8|) 
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is indeterminate) and letting 9 vary from —it to tt. According to fl35|) the wall profile of the 
field 9 as a function of the orthogonal space coordinate x will be given by the relation 

cot {^} =8illh (iL}. (50 ) 

in which the effective wall thickness scale 5* is determined by the relevant mass scale m* as 

^ = — 2 , m* = V if . (51) 

In this case there will be no contribution from the right of (j48j) . so the complete integral 
£7{oo} giving the wall tension T\ works out as 

T I = 4m^ = 2VJ i . (52) 



To get a compound wall trajectory (of the kind illustrated by the curves AC in Fig. H|) 
the variable \ is no longer held fixed but is also allowed, like 9 to vary from — tt to tt. Such a 
trajectory is obtainable by choosing ± = + in equation (1481) . which is then easily integrable 
to give 

cot ( — 1 = k cot | 1 , (53) 



4 J [4 

where k is a dimensionless constant of integration. It takes the value k — 1 in the trivial case 
of the straight diagonal trajectory given by the equation \ = @- F° r an y non-zero value of k, 
the separate energy contributions from the pair of terms on the right of (1491) will integrate to 
the same final result. Thus the various trajectories given by different values of k in fl53l) are 
energetically degenerate, all giving the same compound wall energy 

r I = 4V*5*. (54) 

The condition of being just twice the simple wall energy is an obvious consequence of the fact 
that although their topology is globally interwoven, the separate field combinations 9 and x 
behave locally as a pair of decoupled scalars, so their wall configurations can travel through 
each other without interaction. 

According to the analysis in the preceding sections, this feature 

r n /Ti = 2, (55) 

is just what is marginally needed to ensure stability (against disintegration into pairs of 
Y junctions) of the X-junctions, whatever their crossing angle a may be (not just when 
it is near a right angle as would be needed in models providing a lower value of T^/Tj.) 
Without going into these quantitative details, the stability, as previously claimed [6], of the 
X-junctions in this particular case was heuristically obvious in advance from the lack of any 
mechanism of disintegration (into pairs of Y junctions) in view of the effective absence of 
interaction between the separate 9 and \ fields that respectively characterise the intersecting 
walls. (For sufficiently small but non-zero values of the parameter e in the corresponding 
topologically simple extended model [6] described in the appendix, continuity implies that 
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X-junction stability will still hold except for correspondingly small - meaning highly acute - 
values of the intersection angle a.) 

Although there is no interaction between walls produced by the separate variation of 9 and 
X, it is important to notice that there will be an interaction between the parallel walls in the 
periodic solution produced by the variation of a single one of the separate variables, 9 say, 
for a finite value of the wavelength £ as given by the formula (124"|) for the separation distance. 
In such a lasagne type configuration there will be an interaction, interpretable as an effect of 
mutual repulsion, due to the build up, as £ decreases, of the corresponding orthogonal pressure 
P±. In the separable case under consideration, it can be seen from (|24"1) and (1231) that, in terms 
of the "first" and "second" kinds of elliptic integral [23] , 

K {//} = f 1 . dt E{fi}= /\/ 1 ~ // f dt, (56) 



the relations between pressure, separation, and the enthalpy density H will be given by 

e i — ui \ H E M v. < m 

Wt = V»KW, tTIJu' '^kT2K' (57) 

where 5* is the effective wall thickness as given by (1511) . 

Since the asymptotic behaviour of the elliptic integral is known [23] to be characterised by 



the condition that K{fx} + \n{y/l — /i/4} — > as // — > 1, it can be seen to follow that in the 
long wavelength limit £ — > oo, the relevant pressure will be given asymptotically by 

Pi~8Ke-^, (58) 

while in terms of the tension T\ of an isolated wall, as given by ( 1521) . it can be seen from (125]) 
that the surface energy density per period will be given by the asymptotic formula 

U^T 1 (l + 4e-^) . (59) 

According to (|27|) the velocity of large lengthscale longitudinal perturbations will be given, 
in this long wavelength limit £ 3> <5*, by an expression of the corresponding form 

v± ~ 2 1 , (60) 

a* 

whose reality is what guarantees stability against bunching. 

6. Beta monovac and multivac winding models 

The fields <p and ip in the pentavac model on the preceeding section were set up as the 
respective phases of unimodular complex field variables e 4 * and e 1 ^ that could themselves be 
considered to have been obtained by imposition of a low energy restraint on corresponding 
complex variables $ and ^ whose modulus at much higher energy would no longer be re- 
strained. A new model that would be physically indistinguishable in in the low energy limit 
characterised by the unimodularity condition could be set up by instead taking the inde- 
pendent fields to be the unimodular complex variables e ld and e iX that are defined by the 
combinations 9 and x specified by ( )38l) 
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Whereas the toroidal space occupied by the original pair of unimodular field variables e l< ^ 
and e 1 ^ can be covered by any one of the four large square patches bounded by dotted lines 
in Fig. HI the (much smaller) toroidal space occupied by the new pair of unimodular field 
variables e 10 and e ix can be covered just by the small square with vertices at the positions 
marked A,B,C,D, which in the new model are to be identified, so that there will now be just 
a single vacuum state instead of five. In the new "monovac" model (a doublet generalisation 
of the singlet sine Gordon equation) that is obtained in this way, a simple wall of the kind 
exemplified by the trajectory AB will still be a topologically stable membrane defect, but of the 
kind known specifically as a "winding", to distinguish it from the more commonly discussed 
kind of "open kink" for which (as in the pentavac case) the vacuum states on either side are 
of distinct varieties. 

Both the original pentavac model and the new monovac model can be regarded as special 
separable cases within larger families of respectively pentavac and monovac models for which 
the potential is given in terms of a fixed index (3 > by an expression of the form 

V = ^ (cos 2 {£/2} + cos 2 { X /2}) /3 , (61) 

which includes the special separable example fj40|) as the particular case for which (3 = 1. 

It is apparent that - for sufficiently large values of the angle 5 of deviation from orthoginal- 
ity in Fig. Q]- the possibility of instability of an X type junction with respect to disintegration 
into Y type junctions [9] will occur in such (monovac or pentavac) models if and only if the 
strict inequality /3 < 1 is satisfied. By comparing the double wall trajectory parametrised by 
the symmetric relation \ = $ with the simple wall trajectory given by the fixed value x = 71 
it can be seen from (IBTj) that the energy Tj of the former (as represented by the diagonal AC 
in Fig. |4]) will be related to the energy Ti of the latter (as represented by AB) by 

T n = (v / 2) 1+/3 T I . (62) 

The actual value of the simple wall tension in such a case can be seen to be given in terms of 
the thickness parameter 5* defined by ( |5T|) . and an order of unity factor Ip, by 

ri = 4m*V*, Ip=C^=M=, (63) 

Jo vi — y 

which in the particularly interesting case (3 = 2 gives T\ = nm^r]. 

It evidently follows from (1621 that for /3 > 1 the double wall configuration (such as AC) 
will be unstable with respect to decomposition into a pair of simple walls (such as AB and 
BC) between which there will be an effective repulsion. On the other hand, for (3 < 1 a 
complementary pair of simple walls (such as AB and BC) would be attractive, and in such a 
case, according to (jlj) an X junction will be unstable if but only if the deviation angle 5 in (CQ) 
is large enough to satisfy the strict inequality 

sin 5 > 2 13 - 1 , (64) 

a condition that would evidently be impossible for (3 > 1. 

The foregoing considerations remain true when - to be cosmologically realistic - the mono- 
vac doublet model is extended to a monovac triplet model by the inclusion of a third unimod- 
ular scalar field, e Jcr say, acting in the same way as the others, so that the complete kinetic 
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metric will be given in terms of a relatively large mass scale r/ by 

dA 2 = 7? 2 (d# 2 + dx 2 + da 2 ), (65) 
while in terms of a smaller mass scale m+ the complete potential function will take the form 

V=^ (cos 2 {£/2} + cos 2 { X /2} + cos 2 {a/2})^ , (66) 

which reduces to (I6ip wherever cos{<r/2} = 0. 

It is to be remarked that the potential functions of the same form can be used to charac- 
terise triplet models with not just one but many energetically degenerate vacuum states, by 
supposing that - instead of being obtained directly as phases of unimodular complex num- 
bers - the variables 9 , x > & > are obtained as linear combinations of quantities <fi , ip , £ that 
actually do have this property, meaning that they are phases of corresponding field variables 
e i<f> ^ e iip ^ g i£ _ mos t obvious possibility is to halve the periodicity in each of the three prin- 
ciple directions so as to obtain an octovac model - meaning one with 2 3 (= 8) vacuum states 
- by setting 

9 = 2<\> x = 2^, a = 2£, (67) 
which entails that the kinetic metric (1651) will be expressible as 

dA 2 = V(d0 2 + d^ 2 + d£ 2 ) . (68) 

In the same way, dividing the periods by three instead of two, one can obtain a triplet model 
with 3 3 (= 27) vacuum states by setting 

9 = 3<f> x = 3^, o = 3f, (69) 

for which the kinetic metric (IBT)]) will be expressible as 

dA 2 = 9r] 2 {d^ 2 + d^ 2 + d£ 2 ) . (70) 

An alternative 27-vac triplet model with the same kinetic metric (1701) is obtainable - as a 
triplet analogue of the pentavac doublet model (13"8"j) discussed above - by replacing (1611 by a 
"tilted" (but still conformal) transformation relation of the form 

= 0-2^ + 2^ x = 20-^-2£, a = 20 + 2^ + £ • (71) 

The property of providing same number of distinct vacuum states as the more trivial model 
characterised by the simple transformation (169]) follows from the easily verified condition that 
the Jacobean determinant of the tilted transformation fITT]) has the same value, namely 27. 

By a mutually orthogonal superposition of the lasagne type configuations corresponding 
to variation of 8, \ and a respectively, it can be seen that for (3 > 1 such models can provide a 
robustly stable cubical (rather than merely square) lattice of the kind needed for the averaged 
stress tensor to be isotropic. For a single one of these lasagne configurations, corresponding 
just to the variation of 9 say, with the other two variables fixed at the energy minimising 
values cos 2 {x/2} = cos 2 {<r/2} = 0, it can be seen from (|24|) that for /3 > 1 - instead of falling 



16 



of in the exponential manner described by fl58l) - the dependence on the separation distance I 
of the orthogonal pressure will be given asymptotically, as £ — > oo, by a power law of the form 



m 4 /x \ 2/3/(0-1) 




7 > (72) 



in which 5p is of the same order of magnitude as the original thickness scale 5* defined by 
(15T1) . to which it is related by the specification 

roo 

6p = J P 6 if , J, = J q 7 =====, (73) 

which, in the quartic case (3 = 2 provides the expression J 2 = 4r 2 {5/4}/v^F = in terms 
of a numerical factor close to unity given by <j ~ 1.18. 

It follows that the analogue for (3 > 1 of the asymptotic formula ( |59i) for the surface energy 
density of each slab of thickness i will be given by 

U~T l+ j^\iP± , (74) 

and that according to (1271) the analogue of (160]) for the squared velocity of large lengthscale 
longitudinal perturbations takes the form 

which is manifestly positive, as required for stability, whenever (3 > 1. 

7. Macroscopic comportment 

One of the objections raised against the idea of attributing the cosmological acceleration 
to a regular domain wall lattice was the lack of a plausible mechanism to prevent parallel walls 
from drifting together in the long run and eventually undergoing mutual annihilation. It is 
therefore to be emphasised that this is not a danger in the kind of monovac model advocated 
here, in which the parallel wall number density is interpretable as a conserved topological 
winding number, and the parallel walls are protected against getting too close by the effect 
whereby, when the separation distance I gets too small, the orthogonal pressure P± in (I24|) 
will cease to be negligibly small, and will build up so as to produce a mutually repusive force. 

The simplest cosmologically viable category is that of the separable -(3=1- kind of 
monovac triplet model (effectively a superposition of three independent sine-Gordon singlets) 
which provides a cubic lattice with a cosmologically comoving lengthscale I and hence cell 
number density n = l/£ 3 for which the relations (|5T|) implicitly provide a corresponding 
equation of state whereby the cosmological energy density p , will be given by 

P = 3y, U = H-P ± £, (76) 

while the corresponding isotropic cosmological pressure P , which will be given simply by 

P = P±-2j = 3P ± -2j, (77) 
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which means that the cosmo logically important ratio w = P/p will be given by the formula 



w = - 2 - + ^ . (78) 
3 p 

At a much earlier epoch, when the winding wavelength i would have been comparable with 

or even much less than the wall thickness scale 5+ given by (151]) . it can be seen from (1291) that 
one would have had 

1 . , P ± 1 6ttV , x 

w ~ -- with — ~ - , p — . (79) 

However at the present epoch it is to be presumed that the (cosmologically comoving) length- 
scale £ will have become very much larger, t 3> 5*, and hence that it will be possible to neglect 
the final term in ( 17H1) . which would fall off rapidly with a negative power law dependence on 
£ for (3 > 1, and with an exponential dependence given by the asymptotic formula 

— ' (80) 

for the case j3 = 1, so that (in all these cases) one would ultimately attain the observationally 
admissible [21] value given, using ( 163|) . by 

w~--, with p w 10—^ . (81) 



For a simple perfect fluid such an equation of state would of course be unacceptably 
unstable. The mechanism for the stabilisation confirmed in the present work is of exactly the 
kind that is describable by treating the large scale averaged system as a perfectly elastic solid 
(in which the three relevant fields 9, x, ® would be comoving "base space" coordinates) of the 
kind originally envisaged by Bucher and SpergelpQ, of which it is thus a prefect example. 

In order to have become dominant, the density (I8ip must have recently reached the order of 
magnitude of the baryonic mass density which is given in terms of the proton mass m p and the 
cosmological temperature by pb ~ 10 _8 m p 3 . Since the temperature has a contemporary 
value given in terms of the electron mass m e by C ~ 10~ 9 m c , it can be seen to follow that the 
ratio of the contempory value £ c of the mesh spacing £ to the value (in the millimeter range) 
of the contempory thermal wavelength 0" 1 must be given by 

C 4 ^ 10 27 ^ (^V . (82) 
m p \m e / 

According to the line of reasonning in the preceding analysis [6], the Kibble type wall 
formation mechanism that was envisaged would be likely to provide values of this length ratio 
Q c £ of the order of 10 14 , which would be obtained for m* w 10~ 5 p with p pa m e , or for 
~ 10~ 2 ?7 with p pa 10~ 2 m c . These values are comparatively small (not comparable, as 
was stated due to a transcription error) with respect to the interstellar distance scale, of the 
order of several parsecs, which is what would be obtained for £ c if one used the rather larger 
mass values m* pa p pa 10 _1 m c that were suggested by an analysis of the cruder kind used 
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in earlier work [2]. However it now seems that all such reasoning, whether in its earlier form 
or in the more refined [6] version, should be considered to be effectively obsolete, because 
its generic outcome will be the formation of disorganised lattice configurations that would be 
insufficiently stable, unlike the superposed lasagne type configurations envisaged here. 

The whole question of the wall formation process needs to be entirely reconsidered in the 
context of the stabilisation mechanism considered here, which does not work for random wall 
lattice configurations, but is essentially dependent on the existence of effectively conserved 
winding numbers with sufficiently large values, representing the effect of systematic winding 
in the same sense over a cosmologically large scale. Even for an underlying field model of 
the appropriate kind (as illustrated by the examples considered above) the ordinary kind of 
Kibble mechanism considered so far [21 [6] would give rise to random combinations of positive 
an negative windings, corresponding to what might be called walls and anti-walls, which in 
the long run would be unstable with respect to mutual annihilation. The work in the following 
section suggests, however, that it be possible to overcome this problem by incorporating the 
effect of the kind inflationary process [25J that has been developped for dealing with what 
is known as the horizon problem. On the basis of this revised analysis it will be found 
that the most plausible mean mass values are not smaller but larger than those that were 
originally envisaged, with the implication that the present mesh lengthscale £ c should actually 
be expected to be of intergalactic order. 

8. An inflationary mechanism? 

The horizon problem was posed in the context of traditional cosmology by the need to 
account for the approximate homogeneity - indicative of causal contact in the past - that is 
observed all the way out to the present value of the Hubble radius Ru as defined, in terms of the 
proper time derivative a of a comoving lengthscale a say, by 1/ Ru = a/a. This suggests that the 
radius, Rc say, of effective causal correlation must be substantially larger than this, Rc 3> -Rh- 
The problem is that in the traditional scenario of a decelerating radiation dominated universe 
this Hubble radius has the same order of magnitude as this "particle horizon" radius Rc. To 
get over this, the idea put forward by Guth [25] was that such a restriction would no longer 
apply if, at some stage in the past, there had been a sufficiently long period of inflation during 
which the universe was subject to acceleration, a > of the kind that seems to have recently 
started again. 

When such considerations were originally introduced by Guth [25] , the intended application 
was to the breakdown of grand unification leaving a relic distribution of monopole particles 
whose density per Hubble volume was assumed to have been at least of the order of unity at the 
relevant initial time t h implying an initial number density n- x > R^f. The monopole problem 
was actually even worse than commonly supposed since such a lower limit can plausibly be 
strengthenned by taking account of the consideration that thermal fluctuations corresponding 
to the temperature 0; (in Planck units) at that epoch might have initially produced roughly 
of the order of one of the relevant particles or antiparticles per thermal volume O -3 , of which 
there would have then have been N ~ (-Rci@i) 3 in the correlated volume under consideration. 
As the causal correlation radius increased, the created number in the corresponding volume 
would also have increased by a factor (_Rc/-Ra) 3 ( a / a i)~ 3 to give 



Due to their causal interaction, these particles and antiparticles would however have under- 




(83) 
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gone mutual annihilation, except for a relatively small excess with the same sign (e.g. with no 
antiparticles) whose number would have had a root mean square value that can be estimated 
by a random walk argument to be of the order of \/iV, as given by the estimate 



'N « {RcQ/Zif 2 . (84) 

in which the thermal evolution formula 

0/0 i = Z i a i /a, (85) 

has been used to define a reheating factor Z\ factor of the kind that was supposed by Guth to 
have become extremely large, and that for consistency with the second law of thermodynamics 
must at least increase monotonically from an initial value Z\\ = 1. In terms of the amplification 
factor 

K = R c /Rh, (86) 

this means that the order of magnitude of the number density with which they would have 
ultimately emerged would be expressible as 

/ \ 3/2 

B "(szSfe) ■ (87) 

Since the ratio of the present cosmological radius Rn c ~ 10 59 to the thermal lengthscale given 
by the present cosmological temperature C ~ 10~ 31 is given roughly by _Rhc©c ~ 10 28 , the 
implication is that the number of such particles in a Hubble volume would be given now by 



UcRhc UfPJ ~ 10 ' (88) 

in which the lower limit on the right is obtained by postulating the order of unity values 
that are the smallest conceivable possibilities for the inflation amplitude N c and the reheating 
factor Zi C at the present epoch. Guth's idea [25] for the latter, in the context of the monopoles 
predicted by the Grand unification theory (which was more fashionable then than now) was 
to postulate an enormous value Z lc ~ 10 28 in order to get the right hand side of (IBB]) down to 
near the order of unity. 

At the opposite extreme, in the context of the baryon formation problem, it is to be noted 
that although it is rather large, the maximum value on the right of (1881) is barely the square 
root of the value - nearer 10 76 - of Dirac's cosmological baryon number, thus falling far short 
of what would be needed even to make a single star, so it can not be hoped that a mechanism 
of this kind would suffice to solve the problem of the excess of baryons over antibaryons. 

What such a mechanism might conceivably do, however, would be to solve the analogous 
problem of getting the excess of positive over negative windings that would be needed to 
provide the kind of domain wall lattice discussed above, in which the lattice spacing I would 
correspond to a cell number density n = £~ 3 , so that the formula (1B7|) would provide for the 
the relevant contempory values the analogous estimate 

^HcOc » 10 14 , (89) 



V%z~ ( 



which means that, even for very moderate values of the ratios N c and Z- lc , the contempory 
value £ c of the wall spacing now would at least be large compared with the size of the solar 
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system. To avoid having £ c > Ru c , which would mean that the walls would be entirely inflated 
away (in the manner envisaged by Guth for the monopoles of grand unifcation) it can be seen 
that the inflation factor must satisfy 

y/%Z ic < 10 14 . (90) 



Comparing (|89|) with the matching condition (1821) one sees that it requires that the relevant 
mass ratios should to be related by 

^(^V«io-y^z~ (91) 

On the assumption that the potential mass scale is small compared with the kinetic mass 
scale, meaning m* <C rj, it can thereby be see to follow from (|90|) that the lighter mass scale 
must be subject to a rather severe upper mass limit given by 

m* < 10 2 m e . (92) 



An opposing restriction comes from the astrophysical consideration that, to avoid under- 
mining the experimentally well confirmed scenario of chemical element formation in the tem- 
perature regime G ^ m e , the postulated inflation following the relevant phase transition, with 
G; ^ i], must have been finished before this stage, which evidently imposes the requirement 
that the heavier mass scale should be subject to a lower mass limit given by 

r] > m e . (93) 



9. Conclusions 

To compute the stability of an X type junction in the framework of a particular field 
theoretical model it is not necessary to work with the essentially two dimensional geometry 
that is actually involved: it will suffice to consider the effectively one dimensional problems 
of plane wall configurations of simple and double type, in order to obtain the corresponding 
energy densities Ti and Tj whose ratio is needed. A right angle junction (a = ir/2 in Fig. [1]) 
will always be stable so long as 

7VTi>V2, (94) 

but a larger ratio will be needed for stablity of X junctions at other angles. 

For models of the forced harmonic kind it has been shown that the required tensions are 
obtainable as the field space distances between the relevant vacua, as measured with respect 
to the energy metric given by (|3~4"|) . For the separable models such as the pentavac doublet 
model [6] that was subject to question [9], this geodesic method has been used for the explicit 
evaluation of the relevant tensions, and it has been shown that the stability condition is 
indeed satisfied. Finally attention has been drawn to a cosmologically more promising class of 
monovac triplet models for which the stability condition is guaranteed whenever the relevant 
index {3 exceeds unity. 
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At a macroscopic level, the particular models considered here are well described by the 
formalism of an elastic solid medium, but it is to be observed that, convenient though it 
is [I], this feature is not essential for their stability. The winding stabilisation mechanism 
for the superposed lasagne type configurations envisaged here would work just as well if, 
instead of three (the minimum compatible with an isotropic total stress tensor) a higher 
number of independent scalar fields were invoked. In that case one would obtain more than 
triply superposed lasagne type configurations which would be just as stable, but in which 
there would be too many degrees of freedom for the system to behave as a simple elastic 
(or even hyperelastic [2B]) solid unless coupled with something else that ensured the required 
cohesion. The lack of a plausible cohesion mechanism was one of the weak points in the kinds 
of solid lattice scenario that were originally considered. It is therefor to be emphasised that 
this drawback does not apply to the stability mechanism considered here, for which no such 
cohesion mechanism is needed. 

The drawback that remains is, as remarked in the introduction, the need for the system to 
have been created in a configuration that is wound up coherently in the same sense (without 
string winding defects, and with only "positive" walls as opposed to "anti- walls" ) over a 
sufficiently large cosmological scale. This last drawback is serious, but (particularly in view 
of the likely involvement of an anthropic selection mechanism) it should not be considered to 
be automatically fatal. It should rather be taken as a challenge that is comparable with (and 
perhaps related to) the unavoidable challenges posed by the problem of the parity breaking 
that is associated with the observationally well established preponderance of ordinary matter 
over anti- matter [18], and by the - rather less intractible and perhaps more directly relevant - 
horizon problem, which provides a major incentive for the hypothesis of cosmological inflation, 
and also motivates consideration of effects of multiconnectedness |12j . 

The provisional investigation in the preceeding section suggests that it may indeed be 
possible for an appropriate inflation scenario to provide what is needed. An attractive feature 
of this approach is that while the concomitant restrictions (I90l . (|92l) . (|93|) rule out more exotic 
mass values, they can be satisfied in a very reasonable way by postulating that the relevant 
potential and kinetic mass scales have the most obvious orders of magnitude, namely 

m* pa m e , 77 pa m p . (95) 

This would require that the relevant late inflation stage - initiated with temperature 0i at 
or below the B.e.v. level - should have provided a mean inflation factor y/^ c Zi C ~ 10 11 , in 
which Z ic is the reheating factor and K c is the factor (if any) by which the relevant range of 
causal influence exceeds the Hubble radius i?Hc- The ensuing wall mesh scale would then be 

4 « 10" 3 i? Hc , (96) 

which means that, in a double inflation scenario, a previous such mean inflation factor of 10 3 
(due perhaps just to reheating by 10 6 ) would have sufficed to sweep away any grand unification 
monopoles. It is noteworthy that, in such a scenario, the contempory mesh length £ c would be 
in rather satisfactory agreement with the order of the ten megaparsec scale to which (in the 
context of double inflation) attention has been drawn [27] as the intergalactic threshold for an 
enhancement of the power spectrum of the observed structure. This agreement would also be 
compatible with a more extreme mass parameter ratio and higher inflation temperature, such 
as might be given, with ©i above the T.e.v. level, by m* pa 10~ 2 m c with rj pa 10 4 m p . 
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Appendix A: Derivation from extended field models. 

In a universe with appropriate toroidal topology (so as to admit configurations without 
any string defects) the field models considered above might be postulated to have an essen- 
tially fundamental status. However in might be computationally convenient (and in other 
contexts, such as that of Kibble type mechanisms [IS1 [III 02] , more natural) to consider such 
topologically non-trivial field models to be derived as approximate effective models from ex- 
tended field spaces of higher dimension but with simple geometry, in which the relevant total 
potential V to t contains a sufficiently steeply variable confining contribution, V co { say, having 
a degenerate minimum on a submanifold characterised by the vanishing of a set of functions 
T a say, consitituting a reduced field space with the non trivial geometry characterised by the 
intrinsic metric Gy. This means that to linear order in the deviations T a - provided they are 
chosen locally as extrinsic coordinates of geodesic normal type - the total kinetic metric of 
the extended system will be expressible in the form 

dA t ^ t = (Gij + K ija T a ) d$ l d& + 5 ab dT a dT b , (97) 

where 5 a b is just a Kronecker unit matrix and the coefficients G^ and Kij a are respectively 
interpretable as components of the first and the second fundamental tensor of the submanifold. 
The idea is that form of the total potential is such that it can be decomposed as a sum 

V tot = V coi + V ved (98) 

of a relatively slowly varying residual contribution V rca \ and a much more rapidly varying 
(and usually more highly symmetric) confining contribution given to quadratic order in the 
neigbourhood of the submanifold where T a = 0, by an expression of the form 

Vcof = A- a tT a ? b /8e (99) 

where the variables A~l are the inverse components of a moderate valued positive definite 
matrix A ab , and e is a dimensionless constant that has to be taken to be very small so as to 
get the strongly confining limit. 

The value V of V re d on the submanifold where T a = 0, is all that will be needed for deter- 
mining the effective behaviour of the system whenever the available field energy is insufficient 
to excite significant deviations from this submanifold. Thus, in such a low energy limit, the 
system will be effectively describable just in terms of the reduced model constituted by the 
reduced field space with the metric G^. 

It is of interest to consider the effect on the energy measure (1331 of the deviations from the 
reduced model that will occur when the confining parameter e is not exactly zero, in which 
case, to linear order, the residual potential will be given in terms of a set of coefficients V a by 
an expression of the form 

V red = V + V a T a . (100) 

To linear order in e and in the deviations T a it can be seen from (I97j) that in the extended 
model the total energy measure associated with infinitesimal field variations d$\ dT a , across 
a wall will be given by 

dU t 2 ot = 2V(G ij + K ija T a +G ij V a T a /V+G ij A-lT a T b /8eV)d¥d^ + 2VS ah dT a dT b . (101) 
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It can be seen that the first term will be minimised by taking the deviation components to 
be given by the formula 



T a = -AeVA ab Q b , Q b = (K ijb + G ij V b /V)&'&' , (102) 

in which, as before, the prime is used to denote differentiation with respect to the kinetic 
distance measure dA in the reduced space. For such an energy minimising (and thus equilib- 
rium) trajectory, it follows that - to first order in e - the total rate of variation of the energy 
will be given in terms of its limit value U' in the reduced model by 

t ab 



UU = U'{\ - eVA ab Q a Q b ) , U' = V2V . (103) 

The corresponding integrated total energy over a path through the wall will therefore be given 
to linear order by an expression of the form 

U tot =U-eW, (104) 

with 



U 



JVWdX, W = jA ab Q a Q b vVWdX. (105) 



Appendix B: Examples of extended field models. 

The pentavac doublet model considered in Section 5 was originally obtained [6] in just this 
way by taking the limit e — > in a broken U(l)x U(l) model involving a pair of complex 
fields with a with a flat kinetic metric of the standard 0(4) invariant form 

d\ t 2 ot = 5?7 2 (d$ d$ + d^ dF) , $ = |$|e i0 , $ = |*|e^ , (106) 

and with a potential given as a sum of the form (198j) with a U(l)xU(I) invariant confining 
part 

Kof = ^((|$| 2 -l) 2 + (|*| 2 -l) 2 ) , V t = mf, (107) 
together with a weakly variable residual part given in terms of the combinations fl38|) by 

V red = ^(|$| 2 |*| 2 (cos^ + cosx) + 2/(l-e)) , V* = mi = eV i , (108) 

in which the final constant term has been included to adjust the minimum of the total 
potential to the standard value Vtot = 0. The ratio 



(109) 



will act as a small symmetry breaking parameter. In this case the minimum of the confining 
potential V co { that will be used as the reduced field potential consists of the toroidal locus where 
|$| 2 = \1/| 2 = 1, to whose neighbourhood the field will be closely confined when 3> m+ unless 
the available field energy density is relatively large compared with V*. This neighbourhood 
will include the five vacuum positions where Vtot vanishes, which happens wherever cos 9 = 
cosx = —1 on the nearby toroidal locus where |$| 2 = |\1/| 2 = 1/(1 — e) . 
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The same reduced field model and linearised extension can be more elegantly obtained 
from a variant with a confining potential of the same form (11071) . namely 

Kof = ^((|$| 2 -l) 2 + (|^| 2 -l) 2 ) , (HO) 

but with an adjustment of final term in ( 11081) so as to get a residual contribution of the more 
convenient form 

Vred = ^|$| 2 |^| 2 (cos# + cosx + 2) = ^|$| 2 |*| 2 (cos 2 {#/2} + cos 2 { X /2}) , (111) 

which is such as to ensure that the vacua occur exactly on the submanifold of the reduced 
manifold - namely the locus where |$| 2 = l^l 2 = 1- even for finite values of the dimensionless 
ratio (m*/mfl) 4 . 

To linear order in e and T these variants are effectively equivalent, as they are both 
expressible in the same standard form of the kind used in (197|) and (1 1 1 j) by taking the 
coordinates of the reduced system to be 

$i = y/5r)(f), ^ = y/5r]ip (112) 

and by taking the deviation coordinates to be given by 

T 1 = VSr) (|$| — 1) , T 2 = V5 77 (1*1-1), (113) 

which gives V 1 = V 2 = 2Vj\fhr\. This means that the nonvanishing components of the 
first fundamental tensor on the reduced submanifold (where T 1 = T 2 = 0) will be given 
by Gr n = G 22 = 1, while the nonvanishing components of the corresponding second fun- 
damental tensor will be given by /C U1 = /C 222 = 2/y/b r q. According to (11021) this will give 
Q 1 = 2y/Er) (2(j)' 2 + ip' 2 ), Q 2 = 2a/5?7 (0' 2 + 2 , i// 2 ) . Since the only nonvanishing components of 
the matrix characterising the confining potential will be given by A 11 = A 22 = 5?7 2 /8V* , one 
finally obtains the formula 

A ab Q a Q b = ($ l/4 + $ 2,4 + 4)/2^, (114) 
in which it is to be remarked that the variation rates are subject to the normalisation condition 

$i'2 _|_ $2'2 = i _ 

For a wall of the composite kind obtained by taking \ = 0, one obtains V = V ir cos 2 {9/2} 
so it can be seen from (11051) that the ratio of the corresponding values W^n and U n will be 
given by 

Wz/Uj = 2V i ,A ab Q a Q b /3 , (115) 
whereas for a simple wall the average value of V will only be half as much so one obtains 

Wj/U^V+A^QaQt/S. (116) 

In the latter case as exemplified by the trajectory <p = 2'ip obtained by setting x = (AB in 
Figure 2) the variation rates of the independent angles <p and ip will have a two to one ratio, 
so (|114p will give the value A ab Q a Qb = 11 7/50 V* , from which one finally obtains 

Wx/Ui = 117/150. (117) 
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In the case of a composite wall as exemplified by the trajectory <p = i(>/3 obtained by setting 
X — d (AC in Figure 2) the corresponding ratio will be one to three, which leads to the value 
A ab Q a Qb = 241/100V* , from which one finally obtains 

W^/C^ = 241/150. (118) 

The ratio of the corresponding net wall tensions, T : = U\ — eW\ and T"i — — eW^, can 
thereby be seen to be given to first order by 

Tj 62 s , , 

— = 1 . 119 

2Ti 75 V ; 

It can be seen from this by (J7j) that in such a pentavac model the simple X junctions will be 
stable so long as the intersection occurs with an acute angle a that is not too small to satisfy 
the condition that is given in the small s limit by 



Or 




2.57^, (120) 



but destablisation can be expected to occur for moderately large values of a in cases of the 
weakly confined (and more numerically tractable) kind - to which the work of Avelino et al 
[9] was restricted - for which e is comparable with unity 

Such destabilisation by deconfinement does not occur for the analogous generalisation of 
category of monovac triplet models of the kind considered in Section 6, which can be similarly 
obtained as a low energy limit, for m* <C m$, from broken U(l) xU(l) xU(l) models with a 
flat kinetic metric of the standard 0(6) invariant form 

dX 2 = r] 2 {dOdO + dXdX + dEdE) , (121) 

for 

= \0\e ie , X=\X\e ix , E=\E\e ia , (122) 
with a confining term of the form 

V coi = ((|0| 2 -1) 2 + (IAT-1) 2 + (|^r-l) 2 ) , (123) 



The required form fl66l) of the potential V for the ensuing reduced model is then obtained 
by imposing the unimodular retriction on a weakly variable residual part that is taken to be 
given in terms of a fixed index j3 by 



V red = ^ (|6>| 2 cos 2 {#/2} + |X| 2 cos 2 {x/2} + |r| 2 cos 2 {cr/2}) P , (124) 

so that its dependence on the field moduli will be quartic for /3 — 2, and quadratic in the case 
(3=1. It is to be observed that in the special case (3 = 1 this model will be separable, not just 
in the reduced limit but also (unlike the pentavac case) for finite values of the confinement 
parameter e, so that it is obvious that its X junctions will remain stable even for infinitesimally 
small values of their acute intersection angle a, and it can easily be checked directly that there 
will be a cancellation whereby the relevant analogue of the pentavac limit formula ( 11201) works 
out in this case this case simply to be a c = . 
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An intrinsically equivalent reduced model is obtainable from an algebraically simpler al- 
ternative, albeit with not just one but 8 distinct vacuua, that is specifiable in terms of 6 real 
variables, X- u 3^, for i=l, 2, 3, which combine as a set of 3 complex variables, whose phases 
determine the corresponding angles 6 = 2£i, x = 2(^2, ° — 2C3, in the form 

X i + iy i = Z i = \Zi\e lCi , (125) 

with the potential contributions in ( 1981) given (in a form that is purely quartic for (3 = 2) by 

v** = ^ E(^ 2 + y? - 1) 2 , V red = ^(E^if > (126) 

so that unimodularity in the low energy limit is obtained in the same way as before on the 
assumption that m$ 3> m*. Provided the flat kinetic metric is taken to have the form 

dA 2 = VECd^ + dtf), (127) 

i 

this will give a reduced model that will be physically equivalent to the monovac model 
described above - as can be seen by making the identifications X 1 = cos{#/2}, X 2 = cos{x/2}, 
X 3 = cos{cr/3}, in the unimodular limit \Z ± \ = \Z 2 \ = \Z 3 \ = 1, which leads to expressions of 
exactly the same form (1651) . (1661) as before. By setting = ( t , ip — C 2 , £ = C3 > it can b e seen 
that the reduced model obtained in this way is mathematically equivalent to the octovac model 
characterised by fIBTl) . The octovac triplet model differs in principle from the corresponding 
monovac triplet model in that two kinks not just one are needed to obtain a complete winding. 
However this difference would not be perceptible in practice within the restricted framework 
of the reduced model, as there would be no way of telling whether the periodicity of the phase 
variables should be 2ir or 4n, and therefore no way of knowing whether the relevant vacuua 
were really physically distinct or not. 
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Figure 4: Contours of the potential V, using darkest shading at maxima and bright colours 
to distinguish distinct minima - representing vacuua - labelled A, B, C, D, E, showing the 
approximate location of some alternative paths from A to C, in plot of cf) against ip representing 
a (fourfold) periodic covering of toroidal field configuration space of pentavac doublet model 
with flat kinetic metric Gij cx 5ij and energy metric = 2VGij. The large white dotted 
square contains the single covering range < <fi < 2tt, < ip < 2tt, and the small white 
square contains the range < 9 < 2tt, < x < 27r that would provide a single covering for 
the corresponding monovac doublet model, for which (as would be the shown by a black and 
white printout of this figure) the five kinds of vacuum colour A, B, C, D, E would cease to be 
distinguishable. 
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